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Abstract 

The parametric lattice-point counting problem is as follows: Given an integer matrix A e z mxn , 
compute an explicit formula parameterized by b e U m that determines the number of integer 
points in the polyhedron {x e R" : Ax ^ b}. In the last decade, this counting problem has received 
considerable attention in the literature. Several variants of Barvinok's algorithm have been shown 
to solve this problem in polynomial time if the number n of columns of A is fixed. Central to our 
investigation is the following question: 

Can one also efficiently determine a parameter b such that the number of integer 
points in {xeU n : Ax s= b] is minimized? 

Here, the parameter b can be chosen from a given polyhedron Q Q R m . Our main result is a 
proof that finding such a minimizing parameter is ATP-hard, even in dimension 2 and even if 
the parametrization reflects a translation of a 2-dimensional convex polygon. This result is estab- 
lished via a relationship of this problem to arithmetic progressions and simultaneous Diophan- 
tine approximation. 

On the positive side we show that in dimension 2 there exists a polynomial time algorithm 
for each fixed k that either determines a minimizing translation or asserts that any translation 
contains at most 1 + 1/fc times the minimal number of lattice points. 

1 Introduction 

As many combinatorial optimization problems can be formulated as an integer program, also their 
corresponding counting problems can be formulated as the problem of counting integer points in a 
polytope. Although both problems are hard in general, they can be efficiently solved if the number 
of variables is fixed. This was shown by Lenstra I Len83 1 in the case of integer programming and by 
Barvinok |Bar94| for the integer point counting problem, see also |Bar02, BR07|. 

A parametric polyhedron is a set of the form Pi, = {x e R" : Ax ^ b] for some matrix A e U mxn and 
b e U m . The right-hand-side b is the parameter of iV Barvinok and Pommersheim |BP99| extended 
Barvinok's integer point counting algorithm to the parametric case. They describe an algorithm that 
runs in polynomial time if the dimension is fixed and that computes a quasipolynomialwhose value 
at b equals n Z"|. Since then, several authors described alternative and more efficient algorithms 
to compute this quasipolynomial |KV08, VW07|. Effective implementations of Barvinok's algorithm 
have been provided by De Loera et al. |DLHTY04| and by Koppe andVerdoolaege |KV08 |. Applications 
of parametric integer counting can, for example, be found in compiler optimization |VSB + 07|. Other 
very interesting, though not polynomial-time approaches to the parametric integer counting problem 
can for example be found in IBecOO, BDR021 ILZ55) . 
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In this paper we consider the problem of finding a parameter b e Q such that \Pb n Z n \ is mini- 
mized. More precisely, we consider the following decision problem. 

Integer Point Minimization 
Given: A £ Q rax a rational polyhedron Q c R m , fc £ N 
Decide: 3feeQ: |{x£Z" : Ax =S b}\ sS k 

We remark that if k = and n is fixed, then this problem can be solved in polynomial time with a 
technique of Kannan |Kan92|, see also |ES08|. 

Contributions of this paper 

Our first result is a proof that integer point minimization is 7V\P-complete, even if n = 2, i.e. the para- 
metric polyhedron resides in the Euclidean plane and then even if the parametric polyhedra are the 
translations of some convex polygon along the x-axis. In other words, we show that the following 
problem is iVP-complete. 

Polygon Translation 

Given: A £ Q mx2 and b £ Q ra defining a convex polygon P = {xeU 2 : Ax^b}QU 2 and keN 
Decide: 3A, O^A^l such that P + { ~ A ) = {x + ( ~ A ) : x e P} contains at most k integer points. 

Clearly, this is an instance of the parametric integer counting problem with Q Q U m being the 1- 
dimensional polytope Q = {x e U m : x = b- Xa\, s£ A ^ 1}, with a\ being the first column of A. 

Second, we show that there is a polynomial- time approximation scheme for the optimization version 
of Polygon Translation. More precisely, there exists an algorithm that runs in polynomial time 
for any fixed integer k and either determines a minimizing translation or asserts that any translation 
contains at most 1 + 1/ k times the minimal number of lattice points. 

This result combines techniques from the geometry of numbers with classical techniques from 
discrepancy theory. The discrepancy of a polygon is the absolute value of the difference between the 
number of integer points in the polygon and its area. There is a rich literature bounding this discrep- 
ancy, see, e.g. IHux96l , starting with Gauss' circle problem. Gauss | Gau73l investigated the discrepancy 
of a disk of radius R around 0. He showed that this discrepancy is bounded by 0{R), which implies 
that the number of integer points is @{R 2 ) because the area of the disk is nR 2 . Discrepancy bounds 
also exist for polygons |Hux96|, but they involve the length of the boundary. Instead, we bound the 
discrepancy in terms of the lattice width of the input polygon: the number of integer points in a 
polygon of high lattice width is very close to its area. On the other hand, we adapt a technique of 
Kannan |Kan90| to solve instances with thin polygons exactly. 

2 Polygon translation is NP- complete 

In this section, we provide our main hardness result. First, we discuss the relation of arithmetic pro- 
gressions with the polygon translation problem. 

2. 1 Arithmetic progressions and their pulse functions 

The arithmetic progression defined by the triple {a, k, d) is the set A = {a, a + d, a + 2d, ...,a + kd}. We 
say that a function of the form 
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(£i,yi) (n,yi) 



Figure 1: A quadrilateral for a pulse function p{x) denned by the arithmetic progression {0.2, 0.45, 0.7} 
with e = 0.08. The red line segments have length 0.28, 0.53 and 0.78 (bottom up) and the green line 
segments have length 0.12,0.37 and 0.62 (bottom up) . As the polygon is translated by x to the left, the 
number of integer points inside the polygon is 17 + p{x). 

where e > is a pulse function. The next lemma establishes a relation of pulse functions with the 
polygon translation problem. Figure[T]illustrates the construction with an example. 

Lemma 1. Let p : DS —- {0, 1} be a pulse function for the arithmetic progression A - {a,a + d,...,a + kd} 
with parameter < e ^ j and discontinuities only in (0, 1) . That is, a - £ > and a + kd + e <\. 

Lety\,y2 £ Z such thaty 2 -yi = k. Consider a convex quadrilateral P with vertices {£\,yi), [r\,yi), 
{(2>yi)y and{r2,y2)- Ifwehave 

1. (\<r\ and £2 < r 2 , 

2. {£\} = a + e and {£2} = a + kd + e, 

3. {r\} = a- e and{r 2 } = a + kd- e, and 

4. k\{[£ 2 \-\.£i\) andk\{[r 2 \-[r l \), 

then there exists an M eN so that 

t-l\ 



lit 



[0 



+ P)nZ z \ = M+ p{{t}) for t e [0, 1]. 



Proof. Consider the horizontal slice L,- = P n {(x,y) | y = y\ + i} for ^ i ^ k. Using conditions (2) to 
(4) of the Lemma, one verifies that L,- is a line segment [a,-,/3/] x {yi + i} with {a,} = a + id + £ and 
{j6, } = a + id - £. In other words, as the segment sweeps left, an integer point leaves at all times t with 
{t} = {fit} = a + id- £, and an integer point enters at all times t with {t} = = a+ id+e. Taking into 
account that L, is relatively closed, one has 



1\ 9 if a+ id-£ < {t} < a+ id + £ 



+ Li) nZ<\= Mi + 



, ) 1 otherwise 



for some M,- £ N. In fact, M,- = - \a{] . 

Summing over all L,, ^ i ^ k, and using the fact that the intervals (a + id - £, a + id + £) are 
pairwise disjoint, we obtain the claim of the Lemma. □ 

Next we consider a decision problem involving several pulse functions. 

Arithmetic Progression Meeting 

Given: pulse functions p\ p n (encoded as their parameters k^ J \ £ (;) £ Q) 

Decide: 3xeU: p(x) = E" =1 Pj(x) = 
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Figure 2: The construction of P in the proof of TheoremlU 

We delay the proof of the next theorem to Section l2T2l 

Theorem 2. Arithmetic Progression Meeting is NP-hard. 

We are now in the position to prove that Polygon Translation is AfP-complete. In our reduction 
to arithmetic progression meeting, we restrict ourselves to pulse functions whose discontinuities lie 
in the open interval (0, 1). We can always reduce to this special case using an affine transformation of 
the pulse functions, so that this restriction is without loss of generality. 

Theorem 3. Polygon Translation is NP-hard. 

Proof. Let pi,...,p n : R — » {0, 1} be an instance of Arithmetic Progression Meeting. We will assume 
without loss of generality that all discontinuities of the pj lie in the open interval (0,1). Let p = 
L" =1 Pj- The goal is to construct a convex polygon P such that 

|(f ~ 1 J+P)nZ 2 | = M+p(tt})fora]lfeR 



for some MeN. Then P, v 



, and M form an instance of Polygon Translation which is a Yes- 



instance if and only if p\ , . . . , p n is a Yes-instance for Arithmetic Progression Meeting. 

The idea for the construction of P is straightforward. Lemma [TJ gives us a tool for constructing 
quadrilaterals Pi,...,P n corresponding to the pulse functions p\, .. . , p n , which we then stack verti- 
cally to form the polygon P with 2n + 2 edges. This is illustrated in Figure|2] 

Formally, we define the polygon P as the convex polygon constrained by 2n + 2 lines: the 2n lines 
defined by the left and right edges of the Pj, and the lines through the bottom edge of Pi and the top 
edge of P n . We will argue that, given a proper choice of coordinates for the Pj, the resulting polygon 
P satisfies the following properties: 

1. P has 2n + 2 edges: the bottom edge of P\, the top edge of P n , and n edges each on the left and 
right sides, which are extensions of the left and right edges of the Pj, respectively. 

2. P has 2n + 2 vertices: the two bottom vertices of Pi, the two top vertices of P n , and 2{n - 1) 
vertices which are obtained as the intersection points of lines through the left or right edges of 
adjacent Pj. 

3. Any horizontal translate of P contains exactly the same integer points as the union of the cor- 
responding translates of the Pj. 



The last property is the result that we really need for the reduction. The first two properties merely 
guide us along during the proof. 



Figure 3: Intersection of the lines through the left edges of Pj and Pj+i in the proof of Theorem[3j 



We will use the same notation as in Lemma [TJ but with superscripts indicating which polygon 



Pj we are talking about. We choose y[ J+1 ^ = y£ + 1 for all j 
j = 1 . . . n. We then choose 



1 ... n - 1, and yf = yf + k& for all 



[£\ 



[ef\+3j + 2 



The fractional part of the £^P is chosen to satisfy the conditions of Lemma[TJ Observe that this fixes 
the y. and £ . up to an integer translation of P. 

Let us describe how our choice of the l/^'j establishes the first two properties listed above on the 
left side of P. Observe that the slopes of the left edges of the Pj are strictly decreasing: the slope of 
the left edge of Pj+\ is always strictly less than the slope of the left edge of Pj. 

Furthermore, we claim that the lines through the left edges of Pj and Pj+\ intersect in a point 

(x,y) with y~l < y < y[ ;+1) ; that is, they intersect "between" Pj and Pj+\. To see this, consider the 
point {x' , y') where the line through the left edge of Pj intersects the horizontal line y = y[ J+ ^- By our 
choice of the [£. J , we know that 



x' < gf + 3j + 1< Y£f J +3j + 2= L^i ;+1J J • 

On the other hand, if (x",y") is the intersection of the line through the left edge of Py+n with the 
horizontal line y' = yJ , then 



- i#C/+U| 



x"<£ ( l +l) - 



•3(7 + l)<L^ y+1) J-3j 



-2=[£f\ 



This establishes the claim. The situation is illustrated in Figured 

An analogous choice of coordinates is used for the right edges of the Pj, to obtain negative slopes 
with strictly decreasing absolute value, and well-placed intersection points of the lines through right 
edges of adjacent Pj. 

Together, these properties imply that the boundary of P is indeed an extension of the boundaries 
of the Pj. In particular, every horizontal slice of P with integer vertical coordinates coincides with a 
slice of one of the Pj, and vice versa. This establishes also the third properly listed above. Hence 



r -L> 
,0 



+ p l = E l f f n) +P /I = E(M W + P/(U})) = \f^M L A+p{{t}) 
7=1 V u ) j=l \j=\ ) 



= :M 



for all teU, which completes the reduction. 



□ 



5 



2.2 Simultaneous Diophantine approximation and arithmetic progression meeting 

Dirichlet's theorem is is a classical result in number theory about the approximation of a vector of 
real numbers using rational numbers of equal denominator, see e.g. | HW08 , chapter 1 1] . It states that 
given a\, a n eU and Q e N, there exist q e {!,... ,Q} and p\ p n £Z such that 



a; s£ rr for all / = 1. 

; q qQ lln 



This result is best possible up to constants for the general case. However, we may ask whether a 
better approximation is possible for a specific set of numbers. This motivates the following decision 
problem, which was shown to be ./VP-hard by Lagarias | Lag85 1 . 

Simultaneous Diophantine Approximation 
Given: a\,...,a n eQ, QeN, e>0 

Decide: 3q e {1, .. ., Q} such that | qctj - \qa,j\ | =s e for all j = 1 . . . n 

We use [-J to denote the nearest integer (breaking ties by rounding down). For computational pur- 
poses, the aj must be rational numbers, though their denominators will typically be much larger than 
Q- 

Proof of Theorem^ Let a,\,..., a n e Q, Q e N, e > be an instance of Simultaneous Diophantine Ap- 
proximation. We assume without loss of generality that aj e (0, 1). 

We will define an instance of Arithmetic Prog ression Meeting with pulse functions po, p\, . . . , p n 
as follows. We scale numerators and denominators so that the denominators of the a j and e are all 
equal and we denote their common denominator by D. For every j = 1 ... n, let 



Pjix) 







\x~ij\< ■§:+ 2^ f or some ie{0,l,...,\Qaj\} 



1 else 



The intuition behind this definition is that we would like to have pj (x) = if and only if | xa j - \xa,j J | =S 
£ as in the definition of simultaneous Diophantine approximation. The correction term ^ is needed 
due to the strict inequality required in pulse functions. Furthermore, we define 



PoM 



\x — i\ < for some i e {1,2,.,.,Q} 



2D 

1 else 



It remains to be shown that the original instance of Simultaneous Diophantine Approximation is a 
Yes-instance if and only if p = L" =0 Pj nas a root. 

Suppose q e {1, . . . , Q} satisfies | qaj - \qaj\ | ^ e for all / = 1 . . . n. Then dividing by a j yields 

\qajl £ 

\q -|=£ — 

aj a, 

and hence one obtains p{q) = 0. 

Conversely, suppose that p{q) = 0. Define q := \q\. As paiq) = 0, we have q e {1,...,Q}. Further- 
more, \q-q\<jfi. Let ij e {0 \Qaj\} with \q- ij/aj\ < e/aj + 1/(2D). Then 

\q-ij/aj\ = \q + q- q-ijlctj\ 
^ \q-ijl(Xj\-\q- q\ 
^ \q-ij/aj\-l/{2D) 
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From this, it follows that 

\q- ij/oCj\ < e/aj + l/D 

and since a j £ (0, 1) also that 

\qacj - ij\ <e + l/D. 
Since the denominator of a j and e is D one has 

\qccj - ij \ ^ e 

which shows that q is a solution to simultaneous Diophantine approximation. 



□ 



3 A polynomial time approximation scheme 

In this section, we show the following theorem, which implies a polynomial time approximation 
scheme for the Polygon Translation problem. 

Theorem 4. For every k £ N, there is a polynomial-time algorithm which, given a polygon P and a 
direction v £ Z 2 as input, either computes a translate tv + P containing a minimal number of integer 
points or asserts that every translate of P is a{\ + \ I k) -approximation to the optimal solution. 

The intuition behind this result is that when P is small, we can use integer programming in fixed 
dimension and adapt a technique of Kannan |Kan90| to find an optimum. On the other hand, if P is 
large and contains many lattice points, then only a small fraction of them is close to the boundary, 
and hence the discrepancy relative to the average number of lattice points that one expects based on 
the area of P is small. 

This idea goes back to Gauss' circle problem. Gauss showed that as r grows, the number of integer 
points L(r) in a disk of radius r is asymptotically equal to its area, nr 2 . In fact, he gave a bound on 
the error term \L{r) - nr 2 \ that is linear in r. The heart of the argument lies in counting unit squares 
intersecting the disk and showing that only 0{r) of them lie near the boundary. Observe that if one 
transforms the setting of Gauss' circle problem by a linear map, the disk becomes an ellipse E and Z 2 
becomes a general lattice A. Instead of unit squares we now count fundamental parallelepipeds of A; 
the trick is to use the right parallelepiped. 

The dual of a lattice A is A* = {y £ IR 2 : Vie A: y T x £ Z}. The width w y {K) of a convex body K 
along a dual lattice vector y £ A* \ {0} is defined as 

w v (_K) := max y T x - min v T x. 

" xeK xeK 

The lattice width w{K) of K is the minimum over all choices of y £ A* \ {0}. The lattice width and the 
corresponding dual lattice vector can be computed efficiently in fixed dimension. Note that in the 
linear transformation of Gauss' circle problem, the diameter of the disk becomes the lattice width of 
the ellipse E. In fact, the following theorem implies the discrepancy bound of Gauss. 

Theorem 5. For every lattice A c IR 2 and convex body K c R 2 with lattice width at least k ^ 1, one has 

, N _ voiro , _3_yoim w hereN=\KnA\ 
|JV det(A) 1 55 2k det(A)' wnere iv ~ l A [uv l- 

Proof. Using a linear transformation, we can assume that the Lowner-John ellipsoid |Bal97| of K is a 
unit disk centered at some zeU 2 : 

B{z,l)QKQB{z,2) 
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Figure 4: The set 3 s (K) from the proof of Theorem^ 



Let B = (foi,£>2) be a reduced basis of the dual lattice A*. In other words, ||&i|| s H&2II and its Gram- 
Schmidt orthogonalization satisfies 

b 2 = b^ + pb l ,b^Lb l ,\p\^. (1) 

It is well known that b\ is a shortest non-zero lattive vector. Thus it is the lattice width direction of 
disks. In particular, 4||&i || = w{B{z,2)) ^ w{K) ^ k. Let 

^ = |xei 2 : OS bfx < 1 and s b\x < 1} 

be the fundamental parallelepiped of A associated to B and let £P be its closure. We relate the lattice 
points in K to the area of K by centering one copy of 3? at each lattice point, see Fig.HJ 

1— , . 1— 

0»{K) = {KnA)±-0*= (J x±-SP. 

2 MiTnA 2 



Let R = dK ± \&>. We have 



K\RQ&>{K)QKuR. 



Let 8 > be the radius of i.e. is contained in a disk of radius 8. Since K is a convex body, we can 
estimate 

vol{R\K)^vol{{K + B{0,8))\K) s \dK\8 + nS 2 , 

where \dK\ is the length of the boundary. By Lemma [6] below, we have 8 =s 55 Jjfc- Furthermore, 
|&K"| s \dB[z,2)\ =4n because K is convex and contained inB(z,2). 

471 1 7T 1 3 , 

vol(fl \K)^—- + —-^ — vo\(K), 
3 k 9k k 2k 

where the last inequality follows from vol(^f) ^ vol(B(z, 1)) = n and k ^ 1. It follows that 

3 

voK^UQ) s vo\{K) + vol{R \ K) s vol{K) + —vol(K). 

A lower bound follows from an analogous argument, and combining these inequalities with voliSP^K)) ■ 
N-voliS*) = N- det(A) yields the statement of the Theorem. □ 

Lemma 6. LetB = [b\, b 2 ) be a reduced basis of A* andS? = {xeU 2 : OS bfx < 1 and s Z?Jx < 1} the 
associated fundamental parallelepiped of A. Then the diameter d ofS? is bounded byd^ jj^ . 
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Proof. Using the triangle inequality, we bound d ^ || x|| + 1| y \\ , where x and y are vertices of & adjacent 
to 0. In particular, let x be the vertex satisfying b\x - and b\x - 1. Using notation from the proof of 
Theorem[5] we compute {b^x = b£x + pbfx = 1. Since xLb\ and hence x is parallel to b2, we get 



? 1 4 

|X|| 2 = ——^7 ^ - 



2 II 2 3||<bi|| 2 ' 

The inequality follows from {TJ and || b\ || ^ || b2 II . Similarly, let y be the vertex satisfying b\y - 1 and 
b\y = 0. We compute [b^) 1 y= bly- nbTy= -p, and conclude using Pythagoras' theorem: 

9 1 p 2 1 4 p 2 4 1 



1KM1I 11^2 II 2 ll^lH 2 3 ll fo lH 2 SUfoill 2 ' 

The statement of the lemma follows from 2\/4/3 = 2.309 .... □ 

In the second part of this section, we will show how to find an optimal translate when the lat- 
tice width of P is at most a constant. We extend a technique which was introduced by Kannan for 
parametric integer programming | Kan90 1 . Kannan determines the lattice width direction of the para- 
metric polyhedron P^ as a function of the parameter b, and partitions the parameter space according 
to width direction and according to how the respective lattice hyperplanes interact with the bound- 
ary of the polyhedron. In our case, the lattice width direction is the same for all parameters, since 
we only translate the input polygon. We partition the parameter space only based on interactions of 
the boundary of t v + P with the lattice hyperplanes orthogonal to the width direction of P. Our main 
extension is that we encode counting the number of integer points on lattice slices in an integer pro- 
gram, where Kannan's work only tested for feasibility. Our approach is compatible with partitioning 
the parameter space based on the lattice width direction, and hence the following Lemma can be ex- 
tended to even more general 2-dimensional Integer Point Minimization problems, provided that 
the lattice width is bounded by a constant for all possible parameter values. 

Lemma 7. Given a dual lattice vector y eZ 2 \ {0} such that w y [P) ^ k, the optimal translate ofP in 
direction v e Z 2 can be computed in time 2 0(fclo 8 fc) fo 0(1) , where b is the encoding length of P, v, andy. 

Proof. Using a unimodular transformation if necessary, we can assume without loss of generality that 
y-e\. Let us sketch a simple algorithm to compute the number of integer points in a translate tv + P. 
Let us denote (5 = min xef „ + p e\x the first coordinate of the leftmost point in the translate. Let 

Si = {tv + P)r\{xeU 2 : jq = \/3] + i},i = 0...k 

denote the integral vertical slices of tv + P. Note that some of the Si may be empty. For each slice, we 
can compute the lower end a; and upper end bi and write S; = + [a;, fed. It follows that 

k k 

\{tv + P)nZ 2 \ = £ \Si nZ 2 | = £>d - fail + 1 

i=0 i=0 

We will argue that this algorithm can be encoded into a small number of integer programs that allow 
us to find the optimal t . We start by writing down the minimization of the number of integer points 



9 



based on the a; and b,. 

k 

min £ y, 

1=0 

ji ^ Bt -Ai + l 

bt-K Bt bt 
at ^ Ai < ai + l 
Ai,Bi £ Z 

We can obtain 7 = similarly: 

jS = j8 + tPi 
jeZ, 

where we precompute /3 as the value of ft for t = 0. The only remaining task is to encode the compu- 
tation of the a, and fo, given t and 7. 

Suppose that we know which edge of tv + P the point (7 + i, a,) lies on, and suppose that the 
corresponding edge of P lies on the straight line defined by c T x = d. Then a, is defined uniquely by 
the equation 

c r (7+ i,ai) T - c T tv = d, 

hence we can express a ; - as a linear function in t and 7 and add a corresponding constraint to the 
integer program. 

Unfortunately, the point (7+ i, at) does not always lie on the same edge. Let us separate the trans- 
lation of the polygon into a horizontal and a vertical component, because a vertical translation does 
not affect the incidence between edges and vertical lines. As the polygon is translated horizontally 
the point (7 + i, a,) moves onto a different edge of the polygon when a corresponding vertex of the 
polygon crosses the vertical line X\ = 7+ i. Over all points (7 + /,«,) and (7+ i,b{), such an event 
happens n times - once per vertex - for one unit of horizontal movement. 

Hence we can separate [0, 1) into intervals I\,...,I n with the property that the combinatorics of 
incidences between vertical lines and edges of the polygon are constant for all 7 - ^ within each inter- 
val. This allows us to solve one integer program for each of the intervals, each integer program with 
the added constrained that 7- £ Ij for some j, and appropriate constraints computing the a,- and 
bi as outlined above. Together, these n integer programs cover the entire space of possible values for 
t , and we simply take the best solution found among all of them. Each individual integer program has 
0{k) variables, 2 k + 1 of which are integer variables, and can therefore be solved in time 2 0{klo % k) b 0(l) 
using Kannan's algorithm for integer programming | Kan87 1 . □ 

Proof of Theorem^ We summarize the algorithm as follows: 

1. Compute the lattice width and width direction y £ Z 2 \ {0} of P. 

2. If w (P) ^ 4k, compute an optimal translate using the algorithm of Lemma[71 

3. Otherwise, assert that every translate is a (1 + Ilk)- approximate solution. 

The correctness of the last step follows from Theorem [5) Let A be the area of P and let OPT be the 
number of integer points in an optimal solution. Then for every teU: 

|(^ + P)nZ 2 |^(l + 4-M < (1+ -r) ■ —5- ' OPT < (1 + \)OPT. □ 
8k 8fc 1-Jr k 



(Vi) 
(Vi) 
CVi) 
(Vi) 
(Vi) 
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